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Abstract 

A model independent calculation is given for the nucleon isovector electric 
charge radius which is valid in the limit N c —>oo, — > 0, iV^m^ fixed. This 
expression reduces to that of the Skyrme model in the limit N c m n — > oo. 
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The study of baryons in the large N c and chiral limits is an interesting and subtle problem. 
The interest stems largely from the belief that for many purposes N c = 3 is large enough 
so that an expansion of physical quantities in a 1/N C expansion gives at least a reasonable, 
if crude, description of hadronic properties while at the same time the up and down quark 
masses are light enough so that an expansion in terms of m q is also reasonable. However, 
for some important physical quantities the limits m q — > and N c — > oo do not commute 
|T|Q|. That is, if such a quantity is treated as a function of m q , then the 1/N C expansion is 
not uniformly convergent. This implies that for such quantities expansions in 1/N C and in 
m q cannot be consistent with each other. This issue is significant for both theoretical and 
practical reasons. On the theoretical side it helps make manifest the nature of QCD in these 
important limits. As a practical matter, if the two expansions are inconsistent, at least one 
of the expansions must be invalid. Thus, one may learn about places where one expects 
the expansions to fail and, perhaps, one may be able to correct for the failure of one of the 
expansions. 

This letter will explore the large N c and chiral behavior of isovector electric properties of 
baryons. It has long been suspected that some of these properties are afflicted with problems 
associated with noncommuting limits. The evidence for this has been in the Skyrme model 
H in which Adkins and Nappi M found that the behavior of the isovector charge radius 
of the proton in the limit m\ — * goes as m" 1 while in conventional chiral perturbation 
theory one finds the same quantity goes as Log(m 7r ) times a coefficient proportional to N c 
@||. They speculated that this was due to the lack of commutativity of the two limits. 
Here it will be shown that that this speculation is correct and the underlying nature of this 
behavior will be explained in terms of the role played by the A isobar. 

More specifically, it will be shown: (i) that the leading chiral behavior of the isovector 
charge radius is calculable in large N c chiral perturbation theory — a generalization of usual 
chiral perturbation which is valid in the limit m n —>• , N c —>■ oo , Ncm^ — > fixed. It is 
given by 



|r 2 |N)^ = pj- - Sm tan- 1 \+C 



where 5m = M^ — M^ and C is a constant subleading in either 1/N C , m\ or both, (ii) There 
is a model-independent prediction of the isovector charge radius for any Skyrme model or 
other large N c hedgehog model, in terms of physical observables, which is valid in the limit 
— > , N c — > oo with the large N c limit taken first. It is given by: 

(Nk 2 |N), rf = T |4^ (2) 
16 7r f£ m n 

(iii) This model-independent Skyrme result differs from that obtained in usual chiral per- 
turbation theory with the limit — > taken first, (iv) However, the model-independent 
prediction from all Skyrme-type models is equivalent with the large N c chiral perturbation 
theory prediction provided the large iV c limit is taken prior to the the chiral limit. This 
provides support for the conjecture that the Skyrme model correctly reproduces all model- 
independent predictions of large N c QCD. (v)The 1/N C expansion for the isovector charge 
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radius may have a finite radius of convergence. There is significant evidence that N c = 3 is 
outside this radius of convergence, (vi) Analogous results can be obtained for other isovector 
electromagnetic observables such as the higher moments of the electric form factor of the 
nucleon or the E2 N-A quadrupole transition. 

The A — N mass splitting plays a central role in eqs. ([[]) and (|2]). This is not unexpected. 
As has been noted previously [§], the fact that 

M A - M N ~ 1/N e (3) 

can easily lead to noncommutativity of the large N c and chiral limits. The reason for this 
clear. The leading nonanalytic behavior (in m q or ) for the charge radius in chiral 
perturbation theory is due to small energy denominators coming from nucleon plus one pion 
states. However, in the N c —> oo limit, — — > and there are new states with small 
energy denominators, such as A plus one pion states which can contribute to the nonanalytic 
behavior. 

The Skyrme model || and other hedgehog soliton models such as the chiral-quark me- 
son soliton model ||, the chiral or hybrid bag model 0] or the soliton approach to the 
Nambu-Jona-Lasinio (NJL) model j§| are based on both large N c and approximate chiral 
symmetry. As all of these models behave identically for all of the issues discussed here, the 
phrase "Skyrme" model will be used in this letter to denote any of these hedgehog models. 
Approximate chiral symmetry is explicit in the form of the effective lagrangian. The large 
N c nature is built into the method of calculation. These models are quantum field theories 
which in general are computationally intractable. However, in the large N c limit, the models 
are calculable. The large N c limit is essential in two ways. The first is that in this limit the 
problem reduces to a tractable classical field theory. Unfortunately, the classical solutions 
are hedgehog configurations which break symmetries of the theory; to obtain physical re- 
sults one must project onto states with physical quantum numbers. There is a semiclassical 
projection method based on the separation of collective and intrinsic variables |J. However, 
this technique is valid only in the large N c limit. Thus, calculations in the Skyrme model 
implicitly correspond to working in the N c — > limit — i.e. to leading order in the 1/N C ex- 
pansion. The large N c character of the model is also manifest in that the calculated baryon 
properties will be consistent with Witten's large N c scaling rules for baryons fT2 |. 



It should be stressed that the semiclassical treatment of the problem is in an essential 
ingredient of the Skyrme-type models — unless one specifies a calculational scheme, these 
nonrenormalizable models are not well defined. The fact that the large N c approximation is 
implicit in the semiclassical treatment used in all calculations implies that if some quantity 
has noncommuting large N c and chiral limits, then one expects that semiclassical calculations 
in Skyrme-type models to correspond to taking the N c — > oo first. 

The value of the electric charge radius may be calculated in a straightforward way in 
Skyrme models using semiclassical techniques @. The result depends on the details of the 
model. However, it is easy to see that as m n — > the integrals are dominated by the 
contributions from the tail of the distribution. The pion fields in a hedgehog configuration 
may be written as ir a (x) = 7r(r)f„ where a is the isospin direction. In any Skyrme-type 
model as v — ► oo, the pion field asymptotes to a p-wave Yukawa form with a field strength 
fixed by g wNN : 
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7r(r) -> 7r asympt (r) = ^— (m* + 1/r) (4) 

The isovector charge radius is given by 

(rVi = j / drr 4 y (7r asympt 2 + A(r)) (5) 

where A(r) is a model-dependent function which goes to zero faster than 7r asympt 2 and X is the 
moment of inertia. Inserting the asymptotic form into the integral, using the Goldberger- 
Trieman relation |L3| g^NNf-K = 9a^n (which is true to leading order in the pion mass) and 
the fact that 1/1 = 2/3(Ma — M^) (to leading order in the 1/N expansion) immediately 
yields eq. (H) + corrections which are higher order in either 1/N C or m q . It is worth stressing 
that this result is model independent in the sense that it applies to all Skyrme-type models. 

Of course, this result is incompatible with the known chiral behavior of the charge radius 
in chiral perturbation theory JH|: (r 2 )/ = i = — ^§j2 Log(m 7r /A) , where A is a mass param- 
eter independent of m q . This was originally deduced using dispersive methods £§]. In more 
modern language this result can be obtained from heavy baryon effective chiral lagrangian 
methods fll4l,15 . In such treatments the leading nonanalytic behavior in m 2 is believed to 



be completely reproduced by 1 pion loop feynman graphs. The term proportional to 5g\ 
comes from diagram a. in fig. () while the term proportional to unity comes from diagram 
b. To the extent, that we are interested in the large N c behavior, the term proportional to 
unity may be dropped since qa ~ N c . 

The nonanalytic behavior in m n in diagram a. in this standard chiral perturbation theory 
analysis comes entirely from the nearly vanishing energy denominators associated with the 
7r-N states. However, in the large N c limit of QCD the A becomes nearly degenerate with 
the nucleon; the mass splitting goes as 1/N C . This can be seen directly in the Skyrme 



model [|iO,iT| and can also be deduced in a model-independent way by large N c consistency 



conditions ]TB|]. Thus, in an effective lagrangian treatment of the combined large N c and 



chiral limits, one must treat the A as an explicit low energy degree of freedom. Indeed, in 
the large N c limit there is an entire band of low-lying states with I = J. As will be shown, 
however, for the leading nonanalytic behavior in 1/N C and m q of nucleon properties, only 
the N and A contribute. 

Moreover, in the large N c limit the 7T-N-A coupling is completely fixed. It is sufficient to 
specify the coupling in terms of a heavy baryon effective field theory (HBEFT) |14| in which 
the baryon degrees of freedom are treated nonrelativistically. The HBEFT is the natural 
formulation of chiral pertubation in which chiral power counting rules can be satisfied. The 
effective 7T-N-A coupling in the large N c limit is given by 

£ttNA = -J— (Am',m; -^Km9(m,mj) ^m,ro ; Q%^a + h.C.) (6) 

where 

x ± Nia ±_ (1/2 1 3/2 W 1/2 1 3/2 \ 
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The couplings in eqs. (H) and (0) are model independent; they can be obtained directly from 
the Skyrme model using conventional semiclassical techniques. This can also be obtained 
from large N c consistency relations (T^|. Indeed, the large N c consistency relations can be 
used to show that any correction to these couplings is of relative order N~ 2 . 

The notion that the A should be included as an explicit degree of freedom in chiral 



perturbation theory is not novel. It has been argued on phenomenolgical grounds |L4| that 
the A (or more generally the decuplet in the three flavor case) should be included since 
it is light in the real (N c = 3) world. The formal need to include the A explicitly in the 
combined large N c and chiral limits is completely clear and has been seen unambiguously in 
the nonanalytic behavior of vector- isovector and scalar-isoscalar observables 0. 

The leading order nonanalytic behavior in the combined large N c and chiral limits (with 
N c m n fixed) presumably comes from one pion loops with either nucleons or Deltas in the 
intermediate state. Thus, in addition to the graph in diagram (a.) of fig. () one needs to 
include diagram (a). Note that other states in the I = J band (i.e. states with I = J > 5/2) 
do not contribute to the leading nonanalytic behavior since, by isospin, at least two pions 
are required to connect to the I = J = 5/2 state and more to higher states while the 
leading nonanalytic behavior comes from a single pion loop. Evaluation of this loop is 
straightforward in HBEFT. The contribution to the isovector charge radius from this loop 
using the large N c 7T-N-A vetex from eqs. (El) and (ITI) is 



/, r 2,vAioop _ Wa r , /A) 5 9a 5m , J m^-Sm 

(k l>1=1 ~ 8^ L °g^/ A ) + 4^ Jml - 5m" V 



Note that the Log term precisely cancels the Log from diagram a. and the sum of the two 
immediately gives eq. ([[]). 

It is useful to rewrite eq. ([[]) in terms d = 5m /m n : 

(W\N) I=1 =^- 2 S(d) + C (9) 

where C is subleading and S(d) is a suppression factor associated with the fact that the A 
is not degenerate with the nucleon: 



s(d) = -^^=tan- 1 IJ^-^-) (10) 

This suppression factor s(d) is identical to the one that appears in the leading nonanalytic 
behavior of scalar-isoscalar and vector-isovector as shown in ref. @. 

A few observations are in order about eqs. (|9p and (|T0D. The first is that if one takes 
the large N c limit with m^ held fixed then d goes to zero. Since s(0) = 1, it is easy to see 
that one immediately reproduces the Skyrme model result of eq. (0) in the formal iV c — > oo 
limit. This is significant in that it supports the conjecture that the Skyrme model correctly 
reproduces all of the model- independent predictions of large iV c QCD. A second issue is the 
treatment of eq. ( |T0D when d > 1. This is of immediate consequence to the real world in 
which d ~ 2. The basic point is that the expression must be analytically continued. In 
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principle, there is an ambiguity in this continuation since there is a branch cut at d = 1. 
However, this ambiguity is trivially resolved by the realization that the imaginary part must 
vanish for this observable. Thus for d > 1 



s(d) = - - F =L= tanh- 1 ( J d —^- I (11) 

The formal limit of m n — > is easy to take. One obtains (r 2 ) I=l = — g ^ ^ 2 Log(m 7r / <5m) . 
This disagrees with the correct chiral perturbation theory result by an overall factor of 
5g\/ (5g\ + 1). The disagreement, however, is of order N~ 2 relative to the leading order and 
is due to the neglect of diagram (b.) of fig. (). 

Using eqs. (||) and fllTD and the physical values for g^, 5m, f n and m n yields (r 2 )i=i = 



.72 fm 2 + C where C is subleading. The experimental value is .88 fm 2 . Thus, this simple 
leading order result gives the correct result to within about 20% without including effects 
such as vector dominance. 

It is worth stressing that the pure 1/N C expansion (with m n fixed) for any quantity which 
depends on s(d) is highly problematic. The value of d in the real world is approximately 2.1 
and s(<iphys) ~ -47 . This is quite far from the large N c value of unity. Ideally, one might hope 
that the disagreement between s(0) and s(d p ^ ys ) can be explained in terms of higher order 
terms in the 1/N C expansion which is just a power series in d around d = 0. Unfortunately 
the radius convergence of this series is unity. Thus, the physical value of d ~ 2.1 is beyond 
the radius of convergence. This implies that the inclusion of any finite number of subleading 
terms in the 1/N C expansion will not improve the description and will, in general make 
things worse: denoting the Taylor expansion of s(d) up to n th order by s n (d) one finds, 
for example, s°(<iphys) = 1, s 1 (c?p hys ) m —.34, s 2 (d phys ) « 1.8 and s 3 (d p hy S ) ~ —2.06. The 
obvious lesson from this is that the pure 1/N C is not likely to be convergent for the isovector 
charge radius and is certainly not convergent if it is dominated by pionic tail effects. This 
problem is particularly serious for Skyrme-type models — the only systematic calculational 
scheme currently known for these models is the pure 1/N C expansion. 

It is clear that analogous effects will follow for any isovector electric property, such as 
higher moments of the electric form factor or moments of the N-A transition form factor. 
The Skyrme model will correctly reproduce the leading order large N c and chiral behavior 
provided that the large N c limit is taken first. However, the corrections to this leading 
order 1/N C behavior for physical parameters will be large and the 1/N C expansion will not 
converge. 

The distinction between the leading order large N c results and the physical results is 
particularly important in the case of properties of the A or of the N-A transition. The 
difficulty is that for iV c — > oo the A is stable and quanitities such as the A charge radius 
or quadrupole moment have an unambiguous meaning in terms of measureable quantities. 
However, as N c decreases to the point that the A becomes unstable these quanitities become 
ambiguous in the sense that there is no simple model-independent connection between the 
calculated quanitity and experimental observables. 

An important example of this is the N-A electric quadrupole transition matrix element. 
In the large iV c and chiral limits this matrix element with the large -/V c limit taken first, the 
matrix element corresponds to a well-defined observable directly proportional to the isovector 



6 



charge radius. A quasi-model independent relationship between the two |TT| becomes exact 
in the chiral limit. This gives a transition quadrupole matrix element which is much larger 
than seen in typical quark models. However, as one increases d beyond unity the matrix 
element ceases to have any well-defined experimental signature. Thus, it is hard to use the 
quadrupole transition matrix element to distinguish between various models of the baryon. 
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FIGURES 

Fig. 1: Feynman graphs for leading chiral behavior of (r 2 )/ =1 . The cross represents an 
insertion of the th component of the electromagnetic current, the dotted line is a pion 
propagator, the solid line is a nucleon propagator and the double line is the A propagator. 
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